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1.0  Executive  Summary 


We  have  studied  a  variety  of  fundamental  problems  that  are  important  to  advanced  RF  sources 
based  on  nonlinear  transmission  lines  (NLTL’s).  We  computed  the  Green’s  function  on  a  slow 
wave  structure.  This  is  a  prerequisite  for  the  generation  of  high  power  microwave  from  the 
output  voltage  pulse  of  a  nonlinear  transmission  line.  Many  complimentary  issues  were  studied, 
such  as  electrical  contacts  and  the  effects  of  rough  surfaces  on  the  heating  and  on  the 
performances  of  the  device.  Their  technological  significance  crosses  many  disciplines,  and 
several  existing  and  potential  future  collaborations  with  other  institutions  resulted  from  this 
research. 

The  major  novel  results  are  summarized  below. 

1. We  formulated  the  Green’s  function  in  terms  of  the  eigenmodes  for  a  general  slow  wave 
structure.  It  is  validated  by  a  simulation  code  in  a  specific  example.  [1] 

2.  We  provided  scaling  laws  for  vertical  thin  film  contacts  [4,6],  and  showed  that  edge  heating  is 
not  serious  in  the  thin  film  limit  [7], 

3.  We  provided  scaling  laws  for  horizontal  thin  film  contacts  [6,10,1 1,15],  including  a  3 -terminal 
device  [12],  and  showed  that  edge  heating  is  very  serious  in  the  thin  film  limit  [11,14],  We  also 
found  a  voltage  scale,  which  depends  only  on  materials  property  [13],  that  characterizes  the 
robustness  of  the  material  against  runaway  heating  due  to  electro-thermal  instability. 

4.  We  constructed  a  scaling  law  for  the  limiting  current  for  spherical  and  cylindrical  diode  [18]. 

It  is  the  first  successful  attempt  in  90  years  to  provide  a  simple  analytic  scaling  for  the 
numerically-tabulated  solutions  of  Langmuir  and  Blodgett,  which  are  still  in  use  today. 

5.  We  derived  a  new  scaling  law  for  a  ballistic  diode  [21],  which  is  a  key  component  for  the 
“vacuum  transistor”,  a  high-speed  transistor  envisioned  for  the  future. 

6.  We  identified  the  role  of  the  RF  magnetic  field,  and  the  role  of  RF  electric  field,  in  the 
enhanced  heating  of  rough  surfaces  including  bumps  and  troughs  [24],  at  the  same  time  we 
validated  a  powerful  numerical  code  that  was  developed  elsewhere. 

7.  We  investigated  the  effects  of  surface  irregularities,  and  more  generally  of  random 
manufacturing  errors,  on  the  performance  of  traveling  wave  amplifiers.  [26,27,28] 

Collaborations  were  held  with  the  Air  Force  Research  Laboratory  (Kirtland)  [1],  Naval  Research 
Laboratory  [27],  Caltech  [24],  Singapore  [18],  and  industry  [11,14],  Colleagues  from  the  Air 
Force  Institute  of  Technology,  AFRL  (Wright  Patterson),  UC  San  Diego,  U  of  Maryland,  and 
MIT,  recently  expressed  an  interest  in  our  works,  either  using  it  or  in  collaborating  with  us. 
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2.0  Research  Topics 


2.1  Nonlinear  Transmission  Line  (NLTL) 


High  power  RF  radiation  sources  typically  require  an  electron  beam.  There  are  substantial 
advantages  if  the  radiation  can  be  generated  without  using  an  electron  beam,  such  as  by  using  a 
nonlinear  transmission  line  (NLTL)  to  transform  an  input  pulse  into  high  frequency  components. 
The  key  question  on  this  type  of  beam- free  radiation  source  is  the  conversion  of  the  NLTL’s 
output  voltage  pulses  into  radiation.  The  AC  voltage  pulse  excitation  coupled  to  the  slow  wave 
structure  (SWS)  would  generate  an  electromagnetic  wave  or  set  of  waves  on  the  circuit,  which 
could  then  be  coupled  out  as  radiation.  For  a  general  temporal  excitation,  the  Green’s  function 
for  the  SWS  is  of  fundamental  interest. 

We  have  constructed  the  Green’s  function  for  the  SWS  to  show  how  the  SWS  would  respond  to 
an  impulse  excitation  [1],  The  Green’s  function  includes  all  radial  modes,  and  for  each  radial 
mode,  all  space  harmonics.  We  compare  the  analytic  solution  of  the  frequency  response  on  the 
SWS  with  that  obtained  from  a  particle-in-cell  code,  as  shown  in  Fig.  1.  Favorable  comparison  is 
obtained  when  the  first  few  lower  order  modes  are  resonantly  excited.  This  gives  confidence  in 
the  prediction  of  converting  a  pulse  train  into  radiation  using  a  slow  wave  structure. 


b/L  =  6,  h/L  =  6,  a/L  =  5.8,  //„L  =  ji/2 


a) 


Fig.  1.  Frequency  response  of  the  SWS,  obtained  from  a),  analytical  theory,  and  b),  particle-in-cell 
simulation  using  ICEPIC  [1]. 
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2.2  Bulk  contact 


Because  of  the  surface  roughness  on  a  microscopic  scale,  true  contact  between  two  pieces  of 
conductors  occurs  only  at  the  asperities  (small  protrusions)  of  two  contacting  surfaces.  Current 
flows  only  through  these  asperities,  leading  to  contact  resistance.  Electrical  contact  is  a  very 
important  issue  to  high  power  microwave  (HPM)  sources.  It  is  also  critical  to  thin  film  devices 
and  integrated  circuits,  carbon  nanotube  based  cathodes  and  interconnects,  field  emitters,  and 
wire-array  z-pinches,  etc.  Contact  problems  account  for  about  40%  of  all  electrical/electronic 
failures  [2],  We  therefore  spent  very  significant  efforts  in  the  evaluation  of  contact  resistance. 

For  decades,  the  fundamental  model  of  electrical  contact  remains  that  of  Holm  [3],  which  is 
based  on  simplified  assumptions.  We  have  substantially  generalized  Holm’s  a-spot  model  to 
higher  dimensions,  including  vastly  different  materials  at  the  joints,  in  an  attempt  to  quantify  the 
effects  of  contaminants  on  electrical  contacts  [4,  5,  6,  7].  Simple  and  explicit  scaling  laws  for 
the  electrical  contact  resistance  with  dissimilar  materials  are  constructed.  Sample  calculations  are 
shown  in  Fig.  2.  The  model  assumes  arbitrary  resistivity  ratios  and  aspect  ratios  in  the  current 
channels  and  their  contact  region,  for  both  Cartesian  and  cylindrical  geometries  [4],  The  scaling 
laws  have  been  favorably  tested  in  several  limits,  and  in  sample  calculations  using  a  numerical 
simulation  code. 

If  the  length  of  the  conducting  channel  in  Fig.  2  (a)  is  small  compared  to  its  transverse 
dimension,  the  model  is  extended  to  study  vertical  thin  fdm  contact  [7],  as  shown  in  Fig.  3(a). 
Figure  3(b)  shows  the  general  scaling  laws  for  the  constriction  resistance  with  arbitrary 
resistivity  ratios  and  geometric  aspect  ratios.  In  the  limit  of  small  fdm  thickness,  we  found  that 
current  crowding  in  the  vertical  contact  is  far  less  serious  than  the  current  crowding  in  the 
horizontal  contact  (Section  2.3  below). 
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Fig.  2.  (a)  Generalized  a-spot  contact,  in  either  Cartesian  or  cylindrical  geometries.  Holm's  a-spot 
corresponds  to  the  cylindrical  geometry  with  h  =  0,  a  «  b,  a  <<  c,  (b)  and  (c)  normalized  contact 
resistance  as  a  function  of  b/a,  and  pi/p2,  respectively,  for  cylindrical  geometry,  (d)  total  resistance  of 
a  cylindrical  channel  calculated  from  analytic  theory  (solid  lines)  in  excellent  agreement  with  that 
from  simulation  (symbols). 
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Fig.  3.  (a)  Vertical  thin  film  contact,  (b)  normalized  contact  resistance  as  a  function  of  the 
aspect  ratio  a/h2  in  (a),  with  pl/p2  =  100, 1,  and  0.01  (top  to  bottom).  [7] 


2.3  Thin  film  horizontal  contact 


There  is  another  class  of  thin  film  contacts  where  the  current  flow  in  the  base  is  in  the  horizontal 
direction,  parallel  to  the  base  (as  shown  in  Figs.  4b  and  4c),  instead  of  vertical  to  the  base  (as 
shown  in  Fig.  3a).  Thin-film  contacts  of  this  sort  are  also  ubiquitous  to  modem  electronic 
devices.  Once  more,  contact  resistance  is  a  fundamental  limiting  factor  to  their  electronics 
performance,  especially  for  novel  materials  with  extremely  high  conductivity  used  to  build 
conducting  channels.  How  to  accurately  characterize  contact  resistance  is  an  important  issue. 


We  have  developed  extensive  models  for  thin  film  contact,  including  the  effects  of  dissimilar 
materials  [8-15]  for  both  Cartesian  and  cylindrical  geometries.  For  the  Cartesian  geometry,  we 
have  extended  the  model  to  study  a  general  3 -terminal  structure  [12],  where  the  voltages  at  the 
three  terminals  may  take  arbitrary  values,  as  shown  in  Fig.  4a.  We  solved  the  potential  and 
current  flows  for  such  a  general  structure  by  decomposing  it  into  an  even  solution  (Fig.  4b)  and 
an  odd  solution  (Fig.  4c).  The  normalized  contact  resistance  for  both  cases  is  shown  in  Fig.  5.  It 
is  found  that  the  contact  resistance  is  insensitive  to  the  resistivity  ratio  for  a/h2  <  1,  but  is  rather 
sensitive  to  the  resistivity  ratio  for  a/h2  >  1  where  a  is  the  constriction  size  and  h2  is  film 
thickness.  Typically,  the  minimum  contact  resistance  in  the  even  solution  is  realized  with  a/h  ~  1, 
[8]  for  both  Cartesian  and  cylindrical  cases,  as  shown  in  Fig.  5a.  The  current  crowding  and 
current  partition  are  examined  thoroughly.  It  is  found  that  the  hottest  spot,  i.e.  the  location  of  the 
strongest  current  crowding  and  therefore  highest  temperature  due  to  joule  heating,  occurs  at  the 
constriction  comer  near  the  source  side,  which  is  point  B  in  Fig.  4(c). 


(a) 


(c) 


Fig.  4.  (a)  3-terminal  Cartesian  thin  film  contact  structure  of  dissimilar  materials,  with  arbitrary 
terminal  voltage  at  the  three  terminals,  which  may  be  considered  as  a  superposition  of  an  even 
solution  (b)  and  an  odd  solution  (c).  [12] 
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Fig.  5.  Normalized  contact  resistance  for  (a)  the  even  solution,  and  (b)  the  odd  solution  of  thin  film 
contact  in  Figs.  4b  and  4c,  respectively.  [12] 


We  discovered  an  intrinsic  voltage  scale  that  governs  the  electro -thermal  instability  in 
semiconductors  [13].  The  voltage  Vs  is  deduced  from  the  heat  conduction  equation  as, 

K  = 

where  k  is  the  thermal  conductivity  and  oq’  is  the  rate  of  change  of  the  electrical  conductivity 
with  respect  to  temperature.  It  is  clear  that  Vs  depends  only  on  material  properties  and  is 
independent  of  geometry  and  the  operating  voltage.  Vs  measures  the  intrinsic  tolerance  of  the 
material  to  electro-thermal  instability.  Numerical  values  of  Vs  are  consistent  with  the  well-known 
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properties  of  several  common  materials  [13].  For  example,  the  value  of  Vs  for  SiC  is  2  orders  of 
magnitude  higher  than  that  of  Ge,  Si  and  C,  implying  that  SiC  is  among  the  most  robust 
semiconductor  against  thermal  runaway  (a  well-known  fact). 

Most  recently,  we  extended  our  thin  film  contact  model  to  include  a  resistive  interface  layer  [15], 
which  represents  a  typical  metal-semiconductor  contact,  as  shown  in  Fig.  6a.  We  developed  an 
exact  field  solution  for  the  contact  in  Fig.  6a,  accounting  fully  the  spreading  resistance  caused  by 
the  fringe  fields  near  the  constriction  comer,  which  are  generally  ignored  in  the  typical 
transmission  line  model  (TLM)  in  Fig.  6b.  Our  model  is  compared  with  TLM,  shown  in  Fig.  7a. 
When  the  interface  resistance  rc  is  small,  TLM  becomes  inaccurate.  Figure  7b  shows  the  regime 
when  the  total  contact  resistance  is  dominated  by  the  interface  resistance  <0.5)  or  by  the 

spreading  resistance  ( >  0.5)  due  to  current  crowding. 


(a>  Interface 

resistive  layer  C 


E 


\j) 


+K, 


z\ 

£ 


A  I  /», 

zJ  A  ii  ■ 


B 

O 
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Fig.  6.  (a)  Electrical  contact  including  specific  interfacial  resistivity  pc,  and  (b)  its  transmission  line 
model  (TLM).  In  (a),  an  infinitesimally  thin  resistive  interface  layer  is  sandwiched  between  Regions  I 
and  II.  [15] 
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a/h y 

(b) 

Fig.  7.  (a),  Comparison  of  normalized  contact  resistance  from  our  exact  field  solution  (solid  line)  and 
TLM  (dashed  line),  (b)  Contour  plot  of  RslRTc'”ai ,  which  is  the  ratio  of  spreading  resistance  over  the 
total  contact  resistance.  Symbols  represent  experimental  data  in  the  literature.  [15] 


Issues  on  the  measurement  of  contact  resistance  and  spreading  resistance  are  studied  in  Ref.  [14]. 
The  thin  fdm  contact  models  are  also  extended  to  a  contact  of  arbitrary  shape,  and  to  high 
frequency  AC  conduction  in  a  constricted  bulk  solid  [6,  14]. 
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2.4  Classical  and  ballistic  diode 


(a)  Classical  diode 

For  90  years,  there  has  been  no  scaling  law  for  the  maximum  current  (also  known  as  the  limiting 
current)  in  cylindrical  and  spherical  diodes.  The  limiting  current  density  is  only  numerically 
tabulated  in  the  celebrated  papers  by  Langmuir  and  Blodgett  (L-B)  [16,  17].  We  recently  found  a 
novel  new  scaling  for  the  L-B  law  for  both  cylindrical  and  spherical  diodes,  which  is  expressed  as 
[18] 


-.3/2 


4  \2e  E 

Japp  =~£0  J - j=  +  cylindrical  or  spherical 

9  v  ni  yJD 


where  Ec  is  the  surface  electric  field  on  the  cathode  of  the  vacuum  diode.  This  dependence  is  valid 
whether  Ra/Rc  is  greater  than  or  less  than  unity,  where  Ra  and  Rc  are  respectively  the  anode  and 
cathode  radius.  F  is  the  correction  factor,  for  a  cylindrical  diode  is  given  by  [18] 


Rc) 

UJ 

=  exp 

(s-s0)(7s  +  23X8s-417)' 


143742 


-1 


and  for  a  spherical  diode  by  [18], 


K  ’ 

‘(s-s0X9^-37X4s  +  143)‘ 

UJ 

-exp 

42092 

-1 


where  5  =  ln[ln( l+Rc/Ra)]  and  sq  =  ln[ln(2)].  The  above  scalings  are  accurate  to  within  5%  for  the 
entire  practical  range  3x  10'5  <  Rc/Ra  <  500  that  was  considered  by  Langmuir  and  Blodgett  [16,17], 
as  shown  in  Fig.  8. 
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Spherical  Diode 
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Radius  Ratio  Rc/Ra 


(a) 


Radius  Ratio  RjRc 


(b) 


Fig.  8.  Comparison  of  the  new  scaling  Japp  with  the  Langmuir-Blodgett  law,  for  (a)  cylindrical  diode,  (b) 


spherical  diode.  In  the  calculation,  the  inner  radius  either  Ra  or  Rc  is  1  cm,  the  applied  voltage  is  Vg  =  1 
V.  [18] 


11 


lb)  Ballistic  diode 


In  microelectronics,  a  diode  is  typically  implemented  by  p-n  junctions  with  semiconductors  of 
different  doping.  Rectification  can  also  be  achieved  by  Schottky  contacts  between  metal  and 
semiconductor.  In  the  miniaturization  of  these  diodes,  it  usually  requires  multiple  metallic 
electrodes  of  dissimilar  materials,  which  bring  many  challenges  to  the  nanolithography 
technology.  Parasitic  effects,  such  as  undesired  doping  and  contact  resistance,  may  also  be 
introduced  to  the  junction  by  the  metallic  electrodes.  It  is  thus  important  to  seek  different 
mechanisms  to  give  signal  rectification. 

Current  rectification  based  on  geometric  effects  has  recently  attracted  considerable  attentions 
[19,20],  In  particular,  the  asymmetric  I-V  behaviors  of  a  graphene  based  geometric  diode  have 
also  been  demonstrated  by  Monte  Carlo  simulations  using  the  Drude  model,  as  well  as  by 
experiments  [20].  However,  there  is  still  a  lack  of  analytic  scaling  to  systematically  guide  and 
optimize  the  geometric  design  of  such  devices.  We  analytically  studied  the  current  rectification 
of  a  2D  asymmetric  ballistic  constriction  (Fig.  9)  [21]. 


Fig.  9.  2D  funnel-shaped  ballistic  diode.  Current  entering  the  interface  AB(CD)  is  Iab 
probability  for  charge  carriers  at  AB(CD)  that  can  be  transported  to  CD(AB)  is  fi  (fi).  [21] 


(/cd),  the 


Following  Sharvin  and  Wexler  [22],  we  study  the  ballistic  charge  transport  by  using  the  analogy 
to  dilute  gas  kinetics  (where  the  mean  free  path  of  particles  >  the  size  of  hole  on  the  container 
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wall).  We  found  that  there  is  an  optimal  condition  under  which  the  I-V  asymmetry  is  maximized, 
as  shown  in  Fig.  10.  Our  calculations  are  verified  by  Monte  Carlo  simulations,  including  the 
effects  of  a  general  drift  velocity  [21]. 


Fig.  10.  (a)  The  effective  net  current,  and  (b)  I-V  asymmetry  factor  rj  as  a  function  of  the  angle  a,  for 
various  incident  charge  carriers  velocity  distributions  labeled  by  \g\  [21].  In  (a),  lines  are  for 
analytical  theory,  and  symbols  are  for  Monte  Carlo  simulations.  [21] 


2.5  Enhanced  power  absorption  due  to  surface  rougness 

This  is  a  joint  project  with  Caltech  [24],  encouraged  by  AFOSR  Program  Managers  Drs.  Arje 
Nachman  and  John  Luginsland. 

Interaction  of  electromagnetic  waves  with  rough  surfaces  is  an  important  topic  in  science  and 
engineering:  surface  roughness  is  known  to  cause  excessive  electromagnetic  power  absorption, 
and  it  is  one  of  the  major  limiting  factors  to  the  performance  of  radio-frequency  (RF)  cavities 
and  slow  wave  structures  with  an  impact  on  communication  systems,  particle  accelerators,  and 
material  characterization  at  microwave  frequencies.  The  enhanced  power  loss  due  to  surface  and 
interface  roughness  exerts  significant  effects  on  signal  integrity  of  microelectronic  circuits. 
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We  provided  an  accurate  assessment  of  the  enhanced  electromagnetic  power  absorption  caused 
by  small  local  surface  defects  consisting  of  cylindrical  bumps  or  trenches,  where  we  assume  the 
radius  a  of  the  defect  is  much  smaller  than  the  free  space  wavelength  X  (a  «  X)  so  that  the 
surface  can  be  considered  as  locally  flat  in  the  absence  of  the  defect. 

We  first  considered  problems  whereby  a  semi-cylindrical  bump  lies  on  a  perfectly  conducting 
surface  (PEC)  under  various  electromagnetic  incident  fields,  as  shown  in  Fig.  1 1.  By  using  the 
equivalent  waveguide  approach  [23],  we  analytically  calculated  the  power  absorption  due  to  the 
presence  of  the  bump,  by  decomposing  the  incident  fields  into  three  modes:  (E0,  0),  (0,  H0i),  and 
(0,  H02),  where  the  field  components  are  shown  in  Fig.  11.  Figure  12  shows  the  normalized 
power  absorption  of  the  cylindrical  bump  associated  with  the  incident  electric  field  E0,  and 
incident  magnetic  fields  Hoi  and  H02  [24],  The  analytical  calculation  is  in  excellent  agreement 
with  that  from  the  numerical  code,  high-order  integral  equation  method  (HIEM),  that  has 
recently  been  developed  at  Caltech  [24],  Thus,  our  analytical  theory  provided  code  validation  for 
HIEM. 


Fig.  11.  A  small  semi-cylindrical  bump  on  a  perfectly  conducting  surface  whose  local  wave  electric 
field  is  E0  and  local  wave  magnetic  field  Ho  in  the  absence  of  the  bump.  [24] 
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Fig.  12.  Normalized  power  absorption  from  (a),  E0,  (b),  Hoi,  and  (c)  H02,  as  a  function  of  8/a,  for  various 
A/a,  where  8  is  the  skin  depth  of  the  bump,  and  X  and  a  are  the  incident  wavelength  and  bump  radius, 
respectively.  The  analytical  calculations  (solid  lines)  are  in  excellent  agreement  with  those  from  HIEM 
(squares),  along  with  other  method  (Fourier-Bessel  expansion,  red  dots)  [24]. 
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HIEM  code  was  then  used  to  calculate  the  power  absorption  of  bump/trench  on  a  finitely 
conducting  flat  surface,  which  cannot  be  solved  analytically.  Figure  13  shows  the  power 
absorption  and  the  electric  field  distributions  for  the  bump  and  trench  on  finitely  conducting 
surfaces  [24],  The  incident  field  decomposition  gives  unique,  penetrating  insight  into  how  each 
component  (E0,  H0i,  H02)  contributes  to  the  heating  on  bumps  and  trenches,  which  would  be 
virtually  impossible  in  a  full  scale  numerical  simulation  that  includes  the  total  incident  fields. 


Fig.  13.  Normalized  power  absorption  from  Eo,  Hoi,  and  H02,  as  a  function  of  8/a,  and  the 
corresponding  electric  field  distribution  inside  the  surface  defect  (for  8/a  =  0.1,  /./a  =  1000),  for  both 
bump  and  trench  on  a  finitely  conducting  flat  surface  [24]. 
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2.6  Effects  of  random  manufacturing  error  on  TWT 


Surface  roughness  and  inperfect  fabrication  processes  also  have  detrimental  effects  on  amplifiers.  This 
was  studied  within  this  grant.  The  traveling  wave  tube  (TWT)  is  a  key  element  in  telecommunication 
systems,  satellite -based  transmitters,  military  radar,  electronic  countermeasures,  and  communication  data 
links.  Variations  in  performance  due  to  finite  fabrication  tolerances  in  the  manufacturing  process  can 
lower  the  fraction  of  TWTs  that  meet  specifications  and  drive  up  the  cost  of  manufacturing  [25-28]. 
These  errors  produce  proportionately  larger  perturbations  to  the  circuit  as  the  circuit  size  is  reduced,  and 
also  at  high  frequencies.  Their  effects  on  the  small  signal  gain  and  output  phase  are  shown  in  Fig.  14  [26]. 
A  peculiar  feature,  noted  in  [26]  and  [28]  is  that,  in  the  statistical  evaluation  of  a  large  number  of  samples 
with  random  errors  in  the  circuit  phase  velocity,  a  significant  number  of  these  samples  show  an  output 
gain  that  is  higher  than  the  corresponding  error-free  tube.  We  used  two  very  different  approach  to  show 
that  the  standard  deviation  is  much  larger  than  the  deviation  in  the  mean  from  an  error-free  tube,  thereby 
resolving  this  puzzle  [26,28].  The  analytical  results  on  the  expected  gain  and  phase  compare  favorably 
with  results  from  numerical  integrations  of  the  governing  equation  in  the  absence  of  space  charge,  but  are 
found  to  deviate  from  the  numerical  integrations  with  the  inclusion  of  space  charge  effects  [28]. 

The  effects  of  the  reverse  propagating  circuit  wave  was  also  analyzed  [27].  In  Ref.  [27],  we  found  that 
reflections  may  significantly  increase  the  statistical  effects  on  the  gain  and  the  output  phase.  Moreover, 
gain  ripples  are  introduced  [Fig.  15].  Whether  this  gain  ripple  can  cause  absolute  instability  is  being 
studied. 
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Fig.  14:  (a)  Mean  values  of  the  power  and  (b)  phase  at  the  output  relative  to  the  unperturbed  values 
calculated  by  three  different  methods.  The  standard  deviation  ab  measures  the  variations  in  circuit 
phase  velocity  as  a  result  of  random  errors.  [26] 
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f/GHz 


Fig.  15.  Small  signal  gain  as  a  function  of  frequency  for  several  values  of  standard  deviations  of  the 
circuit  pitch  angle  distributions.  [27] 
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